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Abstract 

In the framework of the Thomas-Fermi model at finite temperature, we show that a 
coohng nondegenerate gas of massive neutrinos wih, at a certain temperature, undergo 
a phase transition in which quasi-degenerate supermassive neutrino stars are formed 
through gravitational collapse. For neutrinos in the mass range of 10to25keV/c2, these 
compact dark objects could mimic the role of supermassive black holes that are usually 
assumed to exist at the centres of galaxies and quasi-stellar objects. Astrophysical 
implications and constraints on the neutrino mass are discussed for this scenario. 

1 Introduction 

One of the most tantalizing puzzles of this Universe is the dark matter problem, whose 
existence is inferred from the observed flat rotation curves in spiral galaxies |2[ , the 
diffuse emission of X-rays in elliptical galaxies and clusters of galaxies, as well as from 
cluster dynamics. The discrepancy between baryon abundances derived from standard 
nucleosynthesis calculations and the observed amount of luminous matter suggests that 
most of the baryonic matter in this Universe is nonluminous, and such an amount of 
dark matter falls suspiciously close to that required by the galactic rotation curves. 
However, although a significant component of the dark matter in galactic halos may 
be baryonic |^], the dominant part of dark matter in this Universe is believed to be 
nonbaryonic, such as to make this Universe nearly critical. Many candidates have been 
proposed |Q| , both baryonic as well as nonbaryonic, to explain the dark matter problem, 
but the issue of the nature of dark matter is still far from being resolved. 

In this paper, we will focus on nonbaryonic dark matter consisting of weakly inter- 
acting heavy fermions. We are particularly interested in fermions with masses between 
10 and 25 keV/c^ which cosmologically fall into the category of warm dark matter. 
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The reason for choosing fermions in this mass range is that these could form supermas- 
sive degenerate objects which may explain, without invoking the black hole hypothesis, 
some of the features observed around the supermassive compact dark objects of masses 
ranging from to W^'^Mq which have been reported at the centers of a number 

of galaxies [5-10] including our own [11-15] and quasi-stellar objects [16-19,21]. 

Interpreting this fermion as a 10 to 25 keV/c^ r neutrino, is neither in conflict 
with particle and nuclear physics nor with astrophysical observations [22|. On the 
contrary, if the conclusion of the LSND collaboration, which claims to have detected 
Ufj, — > i>e flavour oscillations (2^, is confirmed, and the quadratic see-saw mechanism 
involving the up, charm and top quarks |2^, is the correct mechanism for neutrino 
mass generation, the niass may well be between 6 and 32 keV/c^ |2^ which falls 
into the cosmologically forbidden range between 93 h'^ eV/c^ and 4 GeV/c^ with 
0.4 ^ h ^ 1. It is well known that such a quasistable neutrino would pose serious 
problems in cosmology, as it would generate an early matter dominated phase, starting 
perhaps as early as a couple of weeks after the Big-Bang [2C, 22, As a direct 
consequence of the existence of such a heavy neutrino, a critical Universe would have 
reached the current microwave background temperature in less than 1 Gyr, i.e. much 
too early to accommodate the oldest stars in globular clusters, cosmochronology and 
the Hubble expansion age [20|. 

It is conceivable, however, that, in the presence of such heavy neutrinos, the early 
Universe might have evolved quite differently than described in the Standard Model of 
Cosmology. In particular, primordial neutrino stars might have been formed in local 
condensation processes during a gravitational phase transition, shortly after the neu- 
trino matter dominated epoch began. The latent heat produced in such a first-order 
phase transition, apart from reheating the gaseous neutrino phase, might have con- 
tributed towards reheating the radiation background as well. Moreover, the bulk part 
of the heavy neutrinos and antineutrinos would have annihilated into light neutrinos 
and antineutrinos via the in the dense interior of the supermassive neutrino stars 
[22, In this context, it is interesting to note that, even within the (homo- 

geneous) Standard Model of Cosmology, the reason why neutrino masses larger than 
4 GeV/cP are again cosmologically allowed, is precisely because annihilation of the 
heavy neutrinos reduces the matter content of the present Universe. Since reheating 
and annihilation will increase the age of the Universe, i.e. postpone the time when the 



Universe reaches the current microwave background temperature [20, 22 1, it does not 



seem excluded that a quasistable massive neutrino in the mass range between 10 and 
25 keV is compatible with cosmological observations [22, 27]. 

In fact, following ref. |2^, we will show in this paper that degenerate neutrino 
stars may indeed have been formed during a gravitational phase transition in the early 
Universe. While the existence of such a first-order phase transition is well established 
in the framework of the Thomas- Fermi model at finite temperature |2£], the mechanism 
through which the latent heat is released during the phase transition and dissipated 
into observable and perhaps unobservable matter or radiation remains to be identified. 
At this stage, however, it is still not clear whether an efficient dissipation mechanism 
can be found within the minimal extension of the Standard Model of Particle Physics 
or whether new physics is required in the right-handed neutrino sector, in order to 
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avoid supercooling of heavy neutrino matter. Thus, we will henceforth assume that an 
efficient dissipation mechanism exists, within or beyond the Standard Model of Particle 
Physics, in order to make sure that neutrinos and antineutrinos can settle in the state 
of lowest energy in a reasonable period of time. 



2 Cold neutrino stars and astrophysical impli- 
cations 

The ground state of a gravitationally condensed neutrino cloud, with mass below the 



Oppenheimer-Volkoff (OV) limit, is a cold neutrino star |22, 27, 28, Q , in which 
the degeneracy pressure balances the gravitational attraction of the neutrinos. In- 
serting the relativistic equation of state for degenerate neutrino matter into the gen- 
eral relativistic Tolman-Oppenheimer-Volkoff equations of hydrostatic equilibrium, the 
upper bound for the mass of such a neutrino star, the OV limit, turns out to be 
Mov = 0.54195 ml^m-^gy^'"^ [0]. Here m^^ = {hc/Cy/^ = 1.22205 • lO^^ GeV/c^ 
is the Planck mass and gi, the spin degeneracy factor of the neutrinos and antineutri- 
nos, i.e. = 2 for Majorana or = 4: for Dirac neutrinos and antineutrinos. The 
radius of such a compact dark object is Rov = 4.4466i?Qy, where Rqv — ^GMqv I 
is the Schwarzschild radius of the mass Mqv- There is thus little difference between 
a neutrino star at the OV limit and a supermassive black hole of the same mass, a 
few Schwarzschild radii away from the object, in particular, since the last stable orbit 
around a black hole has also a radius of 3i?Qy. For nonrelativistic neutrino stars that 
are well below the OV limit, mass and radius scale as MR^ = 91.869 h^G-''m-^{2/g^ f 



22], as follows directly from the nonrelativistic Lane-Emden equation with polytrope 
index 3/2, subject to the usual boundary conditions at the center and the surface of 
the neutrino star. If we generalize the boundary condition to include a point mass at 
the center, the Lane-Emden equation is able to describe a degenerate neutrino halo 



around a baryonic star |22, 30| as well. 

The idea that some of the galactic nuclei are powered by matter accretion onto su- 
permassive black holes is based on general theoretical arguments and the observation 
of rapid time variability of the emitted radiation which implies relativistic compactness 
of the radiating object. However, so far, there is no compelling evidence that super- 
massive black holes actually do exist, as the spatial resolution of current observations 
is larger than 10^ Schwarzschild radii. The standard routine in investigating the nature 
of the dark mass distribution at the centers of active galaxies, is to observe stellar and 
gas dynamics. However, gas dynamics is usually regarded as less conclusive, since it is 
responsive to nongravitational forces such as magnetic fields. 

Because of its relative proximity, the best place to test the supermassive black hole 
paradigm is perhaps the compact radio source Sgr A*, located at the center of our 
galaxy. Observations of stellar motions near the galactic center |13, 14] and the low 



proper motion (< 20 km sec~^) |3^] of Sgr A* indicate, on the one hand, that Sgr A* is 
a massive M = (2.5 it 0.4) x 10^ Mq object that dominates the gravitational potential 
in the inner (< 0.5 pc) region of the galaxy. On the other hand, observations of stellar 
winds and other gas flows in the vicinity of Sgr A* suggest that the mass accretion rate 
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is about M = 6 X 10~^MQyr^^ [^]. Thus, if Sgr A* is indeed a supermassive black 
hole, as many believe, the luminosity of this object should be more than lO^'^ erg sec~^, 
provided the radiative efficiency is the customary 10%. However, observations indicate 
that the bolometric luminosity of Sgr A* is actually less than 10^'' erg sec~^. This 
discrepancy, known as the so-called 'blackness problem' of the center of our galaxy, 
has been the source of exhaustive debate in the recent past, in which the notion of a 
'supermassive black hole on starvation' has been introduced. A much simpler solution 
to this problem could be to give up the idea that Sgr A* is a black hole. 

Indeed, an alternative model for the mass distribution at the galactic center was 
proposed [^4| recently, in which the customary supermassive black hole is replaced by 
a supermassive neutrino star. It was shown that a degenerate neutrino star of mass 
M = 2.5 X W^Mq, consisting of neutrinos with mass nii, > 12.0 keV/c^ for gi, = 4, or 
nil, > 14.3 keV/c^ for = 2, does not contradict the current observational data. This 
can be seen in Fig.l, where the mass enclosed within a given distance from Sgr A* is 
plotted for a M = 2.5 x 10^ Mq neutrino star embedded in a standard stellar cluster 
with a central density po = 4 x IO^Mq/pc^ and a core radius iZcore = 0.38 pc. We thus 
conclude that for certain values of nii, and gu, this model of Sgr A* is consistent with 
the mass distribution data. 

More recently, also the spectrum emitted by Sgr A* was calculated in the frame- 
work of standard accretion disk theory, assuming that Sgr A* is a neutrino star with 
radius R = 30.3 Id (~ 10^ Schwarzschild radii), corresponding to a neutrino mass of 
nil, = 12.0 keV/c^ for gi, = 4, or nii, = 14.3 keV for g^ = 2 [36|. In such a scenario, the 
accretion disk is contained entirely within the neutrino star. The calculated product of 
luminosity and frequency vLy is presented in Fig. 2. The thick solid line represents the 
case of a neutrino star with M /Mem = 4 x 10~^, whereas thin solid line corresponds to 
M/MEdd = 10""^. Here M^dd = 2.21 • 10"^M yr"^ denotes the Eddington hmit accre- 
tion rate. The short-dashed line describes a calculation with a M = 2.5 x 10^ Mq black 
hole, M/MEdd = 10~^ and an accretion disk extending from 3 to 10^ Schwarzschild 
radii. The long-dashed line represents the case where the accretion rate is artificially 
reduced to M/MEdd = 10"^. As seen in Fig. 2, the neutrino star model reproduces the 
observed spectrum from the radio to the near infrared band very well. Thus, as our 
model fulfils two of the most stringent conditions, i.e. the mass distribution |12, 14] and 
the emitted spectrum ||35|] , we can conclude that the neutrino star scenario is consistent 
with the currently available observational data. 

We now turn to the discussion of the constraints that our model puts on the neutrino 
mass. Assuming that the most massive known compact dark object, which has a mass 
M = (3 lb 1) X 10^ M0 and is located at the center of M87, is a degenerate neutrino star 
at the OV limit, where it is close to being a black hole, the neutrino mass must be in 
the range of 12.5 keV/c^ <m^< 17.7 keV/c^ for g^ = 2, and 10.5 keV/c^ <ni^< 14.9 
keV/c^ for c/^ = 4 (Fig.3). 

If we interpret Sgr A* at the center of our galaxy as a degenerate neutrino star with 
mass M = 2.5 x IQ^Mq and a radius R < 30.3 Id, as derived in the analysis by Genzel 
et al. 14| shown in Fig.l, the neutrino mass must be niy > 14.3 hs\ /(? for g^ = 2, 
and my > 12.0 keV/c^ for g^ = 4 (Fig.3). Of course, for neutrino masses close to the 
minimal mass, the neutrino star will be far away from being a black hole. In fact in 
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our scenario, the much shallower gravitational potential governing the accretion disk 
inside the neutrino star was essential for the success of the calculation of the emitted 
spectrum of Sgr A*. The near black hole situation is obtained in the OV limit with 
nij, = 444 keV/c^ for = 2, or = 374 keV/c^ for g,, = 4. 

Heavy neutrinos can of course also cluster around ordinary stars. Let us therefore 
assume that the sun is surrounded by a degenerate neutrino halo From the 

Pioneer 10 and 11 and Voyager 1 and 2 ranging data we know that the dark mass 
contained within Jupiter's orbit is = (0.12ib0.027) x 10~^ Mq and within Neptune's 
orbit < 3 X IQ-'^M©. Of course the Jupiter data should be taken with caution, 
as Jupiter tends to eject almost any matter within its orbit |3^]. Nevertheless, taking 
Jupiter's data at face value, and interpreting the dark matter as degenerate neutrino 
matter ||2^, the neutrino mass limits are 12.6 keV/c^ < my < 14.2 keV/c^ for g^ = 2 
and 10.6 keV/c^ < < 12.0 keV/c^ for g^, = 4. If the dark matter contained within 
Neptune's orbit is interpreted as neutrino matter |2^], the mass limits are rrii, < 15.6 
keV/c^ for = 2, and m,^ < 13.1 keV/c^ for gi^ = 4. In summary, a mass range of 
14.1 keV/c2 <mu< 15.6 keV/c^ for g^ = 2, and 11.8 keV/c^ < < 13.1 keV/c^ for 
g,^ = 4 seems to be consistent with all reliable data. 

3 Neutrino matter at finite temperature and 
the phase transition 

A gas of massive fermions interacting only gravitationally has interesting thermal prop- 
erties that may have important consequences in astrophysics. The canonical and grand- 
canonical ensembles for such a system have been shown to have a nontrivial thermo- 



dynamical limit [p8| , 3£]. Under certain conditions these systems will undergo a phase 
transition that is accompanied by gravitational collapse pO[ . It is interesting to note 
that this phase transition occurs only in the case of the attractive gravitational in- 
teraction of neutral particles obeying Fermi-Dirac statistics: it neither happens in the 
case of the repulsive Coulomb interaction of charged fermions nor does it in the 
case of the attractive gravitational interaction, when the particles obey Bose-Einstein 
or Maxwell-Boltzmann statistics. To be specific, we henceforth assume that this neu- 
tral fermion is the which is presumably the heaviest neutrino, although this is not 
essential for most of the subsequent discussion. We, moreover, assume that the weak 
interaction can be neglected. 

The purpose of this chapter is to study the formation of a neutrino star during a 
gravitational phase transition. For simplicity, let us assume that the neutrino star will 
be sufficiently below the OV limit, so that we can treat this process nonrelativistically 



p9|| . The gravitational potential energy V{r) satisfies Poisson's equation 

AV = AnGmln^, (1) 

where the number density of the r neutrinos (including antineutrinos) of mass can 
be expressed in terms of the Fermi-Dirac distribution at a finite temperature T as 
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Here /«(??) is the Fermi function 

In[r]) = / ■ (3 

Jo 1 + ^ 

and the chemical potential. It is convenient to introduce the normalized reduced 
potential 

' if^-V), (4) 



tti^GMq 

Mq being the solar mass, and the dimensionless variable x = v/Rq with the scale factor 
^ ''"^ J =2.1377 lyr(il^)''%-/^ ,5) 

Equation (1) then takes the simple form 

l(fv_ _ 3 
X dx'^ 2 



r^/^/. , (6) 



where we have introduced the normalized inverse temperature (3 = Tq/T^ with To = 
muGMoi/kRo. In equation (6) we recover, at zero temperature, the well-known Lane- 



Emden differential equation |30|, 22 1 



The solution of the differential equation (6) requires boundary conditions. We 
assume here that the neutrino gas is enclosed in a spherical cavity of radius R cor- 
responding to xi = R/Rq, in order to prevent the gas from escaping to infinity. We 
further require the total neutrino mass to be Mi,, and we allow for the possibility of a 
pointlike mass Mb at the origin, which could be, e.g., a compact seed of baryonic or 
any other compact matter. v{x) is then related to its derivative at x = xi by 

V [xi) = — [v{xi) , (8 

which, in turn, is related to the chemical potential by /i = kTQv'{xi). v{x) at x = is 
given by the point mass at the origin, i.e. Mb/Mq = v{0). 

Similarly to the case of the Lane-Emden equation, it is easy to show that equation 
(6) has a scaling property: if v(x) is a solution of equation (6) at a temperature T 
and a cavity radius R, then v{x) = A^v{Ax), with A > 0, is also a solution at the 
temperature T = A^T and the cavity radius R = R/A. 

It is important to note that only those solutions that minimize the free energy are 



physical. The free-energy functional is defined as [39|, 

(27r7i) 



2mykT kT kT 



(9) 
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where 

V[nA = -Gml I d\'^^, (10) 
J |r — r I 

and 

W[n,] = --Gml [ d^rd\'^^4^}^^. (11) 
2 J \r — r\ 

The chemical potential in equation (9) varies with density, so that the number of 
neutrinos Ni, = M^/rriu is kept fixed. 

All the relevant thermodynamical quantities, such as number density, pressure, free 
energy, energy and entropy, can be expressed in terms of v/x, i.e. 

nu{x) = r^/^/i (P-) , (12) 

Pu{x) = ^^p-'/'h ((3^) = ^eUx), (13) 
rriyRQ^Ti 2 \ X J 6 

F = ^fiN, + hroRl J d^xn,{x)^^^^^—^ - RI j d^xP,{x), (14) 
E = ^fiN, - hnRl J d^xn,{x)^^^^^^^ + Rlj d^xe^^{x), (15) 

S=^{E-F). (16) 

We now turn to the numerical study of a system of selfgravitating massive neutrinos 
with an arbitrarily chosen total mass M = lOM©, varying the cavity radius R. Owing 
to the scaling properties, the system may be rescaled to any given mass. For definite- 
ness, the mass is chosen as rriy = 17.2 keV/c^ which is about the central value of 
the mass region between 10 and 25 keV/c^ [^], that is interesting for our scenario. In 
Figure 4 we present our results for a gas of neutrinos in a cavity of radius R = 100i?o- 
We find three distinct solutions in the temperature interval T = (0.049 0.311)To; of 
these only two are physical solutions, namely, those for which the free energy assumes 
a minimum. The density distributions corresponding to such two solutions are shown 
in the first plot in Figure 4. 

The solution that can be continuously extended to any temperature above the 
mentioned interval is referred to as 'gas', whereas the solution that continues to exist at 
low temperatures and eventually becomes a degenerate Fermi gas at T = is referred to 
as 'condensate'. In Figure 4 we also plot various extensive thermodynamical quantities 
(per neutrino) as functions of the neutrino temperature. The phase transition takes 
place at a temperature T^, where the free energy of the gas and of the condensate 
become equal. The transition temperature Tj = 0.19442ro is indicated by the dotted 
line in the free-energy plot. The top dashed curve in the same plot corresponds to the 
unphysical solution. At T = the energy and the entropy exhibit a discontinuity, and 
thus there will be a substantial release of latent heat during the phase transition. An 
important and currently still open question is, how and to which type of observable 
or unobservable matter or radiation this latent heat, which can be interpreted as the 
binding energy of the neutrino stars, will be transferred. 
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Figure 1: Various models for the mass enclosed within a distance from Sgr A*: (i) a M = 
2.5 X IO^Mq neutrino star plus stellar cluster, with neutrino masses in the range of 10-25 
keV/c^ and degeneracy factors Qi, — 2 and Qi, — 4, (ii) a M = 2.5 x 10^ Mq black hole plus 
stellar cluster, (iii) a stellar cluster only. Note that the curves for rrii, — 12.013 keV/c^ and 
= 4, and m,y = 14.285 keV/c^ and Qi, — 2, coincide. Data points are taken from [13,14] 
and references therein. 



Figure 2: Comparison of the theoretical and observed spectra of Sgr A*. The thick solid line 
corresponds to the case of a neutrino star of mass M = 2.5x IO^Mq with M / MEdd = 4x 10~^, 
while the thin solid line represents M / MEdd = 10^"^. The short-dashed line describes a 
calculation for a M = 2.5 x lO^M© black hole, with M / MEdd — 10"^ and an accretion disk 
extending from 3 to 10^ Schwarzschild radii. The long-dashed hue corresponds to the case 
when M / MEdd is artificially reduced to 10~^. Data points in the < 40 GHz region are upper 
bounds [35] . Note, that the thick solid line fits the most reliable data points with error bars. 



Figure 3: Neutrino mass limits based on various assumptions on the supermassive compact 
dark objects at the center of M87 and our galaxy. Also shown are mass limits derived from 
a possible degenerate neutrino halo around the sun. The dotted u^- mass range is derived 
from the LSND experiment and the quadratic see-saw mechanism. 



Figure 4: Density distribution normalized to unity for condensate-like (solid line) and gas- 
like (dashed line) solutions at T = T^. Free energy, energy and entropy per particle are 
shown as a function of temperature. Temperature, energy, and free energy are in units of 
To. 
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